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In one study (Grosfils et al, 1999) of the propagation in random 
Delaunay lattice particle on arriving at a site in the lattice deflects over 
the largest possible angle to either the right or the left, depending on 
the right- or the left scattering nature of the site. After the particle has 
passed through the site the latter goes into the reverse state. When 
a similar study is done on the triangular lattice, the entire trajectory 
quickly becomes confined to a particular strip which is bounded by 
two adjacent parallel lines of the lattice. They explain the propagation 
as being due to a blocking mechanism which prevents the particle from 
moving in a direction opposite to the propagation direction for more 
than a few steps. 


The flow of fluid in porous media follows the Darcy’s Law which 
states that the rate of flow through such a medium is proportional to 
the potential energy gradient within that fluid. The constant of propor- 
tionality is the hydraulic conductivity, which is a property of both the 
porous medium and the fluid moving through it, v = Q/A = —kdh/dl. 


The average velocity over the entire cross section is called the su- 
perficial velocity, v; = Q@/A = t/pA. The velocity that is based on 
the actual open space within the porous media is called the interstitial 
velocity, vu, = Q/(eA) = tm/epA, where ¢ is the porosity, ¢ = V,/V, = 
(Vi — Vs)/Vi, Vo, Vs and VY are respectively the void, solid and total 
volume. 


The definition of porous media is different from that of porous ma- 
terials. A porous medium is a medium through which other substance 
may pass, whereas a porous material merely means some certain kind 
of material the internal structure of which is filled with pores. Most of 
porous materials can be used as porous media. All of them are useful 
because of their internal structure, and though we do not need to know 
this to be able to use them, we do have to understand it if we want 
to use them efficiently, or if we want to improve upon some of their 
particular properties. 


The outer bark of Quercus suber L., for example, has a property 
that is ideal for its commercial use as cork. Its material property is such 
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that it neither shrink nor expand when stretched or compressed. This 
makes it an excellent insulator of both heat and sound as well as good 
for damping vibration (cf Ashby, 1990). It has high frictional coefficient, 
and is both impervious to liquids and chemically stable. These proper- 
ties make it ideal as wine corks. By understanding its internal structure 
and mechanism that makes it stay in the same shape when experiencing 
external forces, scientists have succeeded in making a synthetic mate- 
rial that shrink, instead of expand radially when pressed. When such 
material is used as cork, it can be easily pushed inside the neck of a 
bottle. When we release the force it will expand to its full size and fit 
snuggly where we placed it. 


Pores introduced into a fabric can also increase its commercial 
value by giving it a lighter weight, higher heat retention rate, which 
makes it ideal for both ski and snow-board wearer and casual wear 
(JETRO, 2001). 


Some filters have a fibrous structure, for example the dust filter 
studied by Bahners and Schollmeyer (1986). 


The filtration coefficient X in —dc/dl = Ac is not a constant but 
change with time because of particles adsorbed by the bed. If one 
assumes that particles are bound to the wall by London-van der Waals 
force only, then A = K(€9 — 0) where o is the specific deposit of solids 
in filter bed described as volume of solid per unit filter volume, €o the 
initial porosity of the bed and K a function of London-van der Waals 
constant, d, «, 7 and W. The equation of continuity is Wo(dc/dl) = 
d(a)/d@. Cake is stabilised by the flow force and consolidates when 
it has reached a critical thickness and the velocity dropped below a 
critical value. The cake pressure is highest at the interface with the 
filter, so it is here that it starts to consolidate. Vibration is normally 
used to loosen it. 


Frederic Lewis (1928) in his study of cucumber cells finds that the 
number of sides in each of them are usually 6. These tall epidermal cells 
cover in a single row the rind of cucumber in a very regular manner 
with the average number of sides touching other cells on the same layer 
is 5.996. When these cells divide, the plane of division usually crosses 
a short axis of the polygon, which implies that the dividing partition 
extends from one wall of the cell to another. The number of sides 
and the number of vertices being equal, if we let no be the number 
of vertices in the original cell while m1 and nz that of the two cells 
(c, and c) thus resulted, then, because each of the new cells has two 
new vertices, n1 + n2 = no + 4. Dividing cells have an average number 
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of sides 6.992. Division of a pentagonal cell invariably produces a 
pentagon and a quadrilateral. Sixty-eight per cent of the hexagonal 
cells division is into a pair of pentagons, but they could also be divided 
into a quadrilateral and a hexagon. A heptagon normally divides into 
a pentagon and a hexagon, but some five per cent of the division gives 
a heptagon and a quadrilateral. Octagons usually divide into a pair 
of hexagons, but 24 per cent of the division creates a pentagon and 
a heptagon, and it is theoretically possible to have a quadrilteral and 
an octago. All the nonagons division observed were into hepagons 
and hexagons. One layer of regular hexagons superimposed on top of 
another give a projection of tetrakaidecahedra. The underlying cells 
become facetted in such a manner that they may become a part of a 
heptakaidecahedron. 


There are similar patterns of grains in a polycrystalline ceramic, 
magnesium oxide, cadmium, etc. The moments of distribution above 
n= 6is pm = >, (n—6)%f, where f, is the fraction of cells with n 
sides. The second moment 2 > 0 unless all cells have six sides, and we 
have a purely topological relation }°,, ninefn = W2 + 36. If n& x 1/ne, 
then n° = (6—a + bu?/6) + (6a 4+ (1 — b)y2)/n. If a = 1 and b = 0, this 
equation is reduced to n§ = 5+ (6+ p2)/n and furthermore if p2 = 0, 
n€ =5+8ne. Ina polycrystal the distribution of ne does not usually 
vary as the grains grow. Typically w2 = 2.4. Soap foams resemble a 
polycrystal (Aboav, 1980), and nf = A+ B/ne, n& = (6—a)+ (6a4+ p2)/n, 
a= 1.2. 


The grains of polycrystals are not arranged at random but in a 
characteristic way which can be expressed in simple terms and seems to 
be scale-free (Aboav, 1970). The average number of sides of neighbours 
of a grain is n° =5+8/n, where n, is the number of sides of the grain. 
Aboav (1970) studies grains of polycrystalline magnesium oxide and 
finds ne = 5.85. He also finds that most of the time n& > n. which 
seems like a contradiction but is because the probability that a point 
lies in a grain of a particular shape depends on the size and abundance 
of such grains. He also finds the average number of n-sided grain 
n2—© x (n—2). Grains of a polycrystal are different from cells of a soap 
foam in that they possess a stable grain diameter, which only depends 
on the temperature. 


This stable diameter of grain, de, is the average grain diameter at 


which the growth ceases, and is de” = c(T — Ty) (Aboav, 1971). For 
cadmium Tp = —53°C, 0°C < T < 170°C. The distribution of grain size 


iS: 2) = 2p, EXP {-o? [(@/en)? = )"} where z is the number of grain 
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sections in a plane section, a a constant, x the diameter of grain section 
and z,, the value of x at zm. 


The shape of individual grains in polycrystal and the inter gran- 
ular surface can be found by stereoscopic microradiography, but using 
random plane section is more convenient in practice (Aboav and Lang- 
don, 1969). 


Gas turns into liquid by the growth of larger and larger clusters. 
But unlike percolation in networks, clusters in condensation process 
are not well defined. Crystalisation is also characterised by the growth 
of one phase within another. But here the orientations are an inherent 
part of the clusters themselves, and there is a long range coordination 
among clusters. One important thing in percolation is for the system 
size to be infinite. Random discs overlapping one another in a contin- 
uum helps correct counts of bacteria cultures (cf Essam, 1980). Let p- 
be the critical probability and P(p) the percolation probability. Then as 
p approaches p, from above, P(p) = (p—p-)® where 0.4 < 8 < 0.5 is 
the critical exponent. If p approaches p, from below, S(p) & (pe — p)~7 
and &(p) © (pe —p)~” where S(p) is the conditional average (s|F), s the 
number of particles in a cluster, and F the event that the origin O is oc- 
cupied and belongs to a finite cluster. In other words, S(p) is the mean 
size of cluster at the origin given that F occurs. This last event occurs 
with probability pr = p(1 — P(p)) and pr = p for p < pe. Furthermore, 
S(p) = pp, C(r,p) and €(p) = [prS(p)|"! 5, r°C(r, p) where the pair- 
connectedness function is C(r, p) = (n(r)|F)pr, and n(r) is the indicator 
defined to be one if r is connected to O and zero otherwise. When F oc- 
curs, s = }>,.n(r). There are two different definitions of the critical prob- 
ability, pe = sup{p|P(p) = 0} and a, = sup{p|P(p) = OandS(p) < co}, 
sometimes denoted by py and pr for Hammersley and Temperley re- 
spectively. By definition, 7, < p,. When p approaches p, from above, 
y' and v' are similarly defined respectively by S(p) ® (p — p-)” and 
E(p) © (p— pe)” It is generally assumed as 7! = y and v'! = v. No 
proofs for these assumptions exist, even though they are consistent 
with the series expansions. Estimates of y and vy are 1.6 < y < 1.7 
and 0.8 < v < 0.9. In a dilute ferromagnet a cluster containing s 
spins each of which has a unit magnetic moment has a probability 
p = exp(5sh)/ [exp(—3sh) + exp($sh)| of being parallel to the magnetic 
field H > 0. Here h = 2H/kpT where kpT is the thermal energy. For an 
infinite cluster, p = 1. The zero-field magnetic moment is puo(p) « P(p). 
The zero-field magnetic moment is po(p) « P(p) ~ (p—pc-)®. The field 
dependent magnetic moment at p, is ((1 — exp(—sh))/(1 + exp(—sh)))r 
and 3(1—G(pe,h)) < fe(h) < 1 — G(pe, h) where G(p,h) = (exp(—sh))x- 
Therefore the critical probability 6 is defined by p-(h) ~1—G(pe, h)hi/?. 
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The correlation function between o on site i and j is defined by Ty; = 
(oi0;)r — {o%)r(o;)r where (-)r is an average over spin states. Then 
Ij; = 1 when i and j belong to the same cluster and zero otherwise. 
The fluctuation formula is kgT(x) = 30 ,(Tij) where (Vij) = C(rj — ri, P), 
the mean susceptibility is (vy) = prS(p)/kseT ~ (p- — p)~7 and the mean 
free energy is (F') = (kgT 1n2)K(p) where K(p) is the mean number of 
clusters and K(p) ~ (p- — p)*~%, where the index assigned to the third 
derivative divergence is 1+ a, a  —0.5. 


Shape is the most fundamental geometrical property. Shape and 
size are all the geometrical information that remain when location and 
rotational effects are filtered out from an object (Dryden and Mardia, 
2002), and between these two you can take the size away so that only 
shape remains for further analysis. When we talk about particle sizes 
in simulation, it is usually the case that we have already assumed 
some kind of particle shape. This is because the definition of size is 
only meaningful if you have some idea about the shape. Shape anal- 
ysis works with landmarks, which are also known as anchor-, control- 
, design-, key-, model-, profile-, or sampling points, facets, markers, 
nodes, sites, fiducial markers, etc. Dryden and Mardia (ibid.) work with 
three types of landmarks, viz. anatomical-, mathematical- and pseudo 
landmarks. Their work could become very interesting if combined with 
another problem of object location (cf Tiyapan, 1996). A landmark can 
be unlabelled or labelled with a name or number. A particular member 
of the shape set which is used as a representative of all other members 
is the icon of that set. For the shape analysis in two dimensions the 
thin-plate spline is a convenient tool which is bijective and is analo- 
gous to the monotone cubic spline. 
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